Long timescale dynamics of spin textures in a degenerate F=l ^' Rb spinor Bose gas 
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We investigate the long-term dynamics of spin textures prepared by cooling unmagnetized spinor 
gases of F = 1 ^^Rb to quantum degeneracy, observing domain coarsening and a strong dependence 
of the equilibration dynamics on the quadratic Zeeman shift q. For small values of \q\, the textures 
arrive at a configuration independent of the initial spin-state composition, characterized by large 
length-scale spin domains, and the establishment of easy-axis (negative q) or easy-plane (positive 
q) magnetic anisotropy. For larger |^|, equilibration is delayed as the spin-state composition of the 
degenerate spinor gas remains close to its initial value. These observations support the mean- field 
equilibrium phase diagram predicted for a ferromagnetic spinor Bose-Einstein condensate, but also 
illustrate that equilibration is achieved under a narrow range of experimental settings, making the 
F = 1 ^^Rb gas more suitable for studies of nonequilibrium quantum dynamics. 



Extensive efforts are underway to utilize ultracold 
atoms as simulators of condensed-matter systems, mak- 
ing use of their broad tunability and amenity to novel 
experimental probes [1 . The success of these endeavors 
hinges on the cold-atom system's arriving at equilibrium, 
matching the conditions of solid state systems that reside 
at or near the thermal equilibrium state. However, given 
their limited lifetime, weak interactions, long dynami- 
cal timescales, and several nearly conserved quantities, 
whether equilibrium is reached in cold-atomic materials 
has become a pressing scientific question [2H4] . 

One promising system for studying materials physics 
is the spinor Bose gas, comprised of atoms that may ex- 
plore all Zeeman states within a spin F manifold. At 
low temperature, quantum degenerate spinor gases ex- 
hibit spontaneous symmetry breaking O [6], magnetic 
order, and intricate spin textures [T, "8^ governed by the 
interplay of magnetism and superfluidity. Several works 
have examined the ground state of a spatially uniform 
F = 1 spinor gas for which the spin-dependent con- 
tact interaction, characterized by a coupling strength 
C2 = 4:7rh'^Aa/m < 0, favors magnetized spin states [9]- 
[12]; here Aa is a difference in 5- wave scattering lengths, 
and m is the atomic mass. Under the application of a 
quadratic Zeeman energy qF^^ with F^ being the dimen- 
sionless spin projection along an experimentally selected 
axis f , and neglecting linear Zeeman shifts, both mean- 
field and exact many-body treatments predict a symme- 
try breaking phase transition between paramagnetic and 
ferromagnetic phases at a positive quadratic Zeeman shift 
Q = Qo = 2|c2|n where n is the gas density. Within the 
ferromagnetic phase, the magnetization is favored to lie 
transverse to x for g > (easy-plane) and along x for 
q < (easy-axis). 

This predicted phase diagram (for q > 0) was sup- 
ported by studies of F = 1 spinor condensates of ^^Rb 
in the single spatial-mode regime, for which the conden- 
sate dimensions in all directions were comparable to the 
spin healing length, defined by ^s = (87rAan)~^/^ [13]- 
[T5] . However, recent experiments on spatially extended 



spinor gases revealed variegated magnetization textures 
produced upon cooling unmagnetized spinor gases into 
the quantum degenerate regime [8 . Theoretical works 
examined whether such textures could be favored at 
equilibrium by the substantial magnetic dipolar interac- 
tions in the ^^Rb gas; however, while some inhomoge- 
neous textures do appear energetically favored, or at least 
metastable, they are predicted to evince spin-density 
modulations only on much longer lengths scales than ob- 
served experimentally p!6H2T] . 

To resolve this discrepancy, we examine quantum- 
degenerate spinor Bose gases over evolution times of 
several seconds, much longer than previously [8]. We 
find that cooling non-degenerate gases with different ini- 
tial spin compositions produces spin textures that differ 
strongly at short equilibration times, consistent with Ref. 
[8 . However, at least at small values of \q\ ^ h x 10 Hz 
and at the longest equilibration times, these textures 
do tend toward a common state characterized by large 
commonly magnetized regions and easy-plane/easy-axis 
magnetic anisotropy consistent with the predicted equi- 
librium phases of ferromagnetic spinor condensates with 
contact interactions. For larger values of |g|, the re- 
distribution of Zeeman populations needed to yield the 
expected equilibrium states is dramatically slowed, pre- 
venting the appearance of a common state from differing 
initial populations. 

For our experiments, we produce optically trapped, 
spin-polarized atomic gases in a manner similar to Ref. 
[22] , with initial temperatures of 30 jaK and initial trap 
frequencies of {oOx^^y^^z) = 27r(220, 4150, 63) Hz in a fo- 
cused, linearly polarized light beam with a wavelength 
of 1064 nm. We next prepare unmagnetized spin mix- 
tures with fractional populations in the three Zeeman 
sublevels, (rnj with mp = {+1,0,-1}, by the applica- 
tion of resonant 7r/2 rf pulses in the presence of a mag- 
netic field gradient [F . The spin mixtures are then cooled 
by lowering the optical trap depth over 2.4 s to a fi- 
nal trap depth of 1.5 /iK and final trap frequencies of 
(ujx^oJy^oJz) = 27r(25,480, 7.3) Hz. We then allow for an 



equilibration period of up to t = 4 s, where t = denotes 
the end of the ramp-down of the optical trap depth. 

Throughout the evaporation and equilibration periods, 
a static magnetic field of B = 267 mG was applied along 
the X axis, with a spatial inhomogeneity of less than 1 
fiG across the extent of the gas. This field causes the 
transverse magnetization of the sample to precess with 
respect to the fixed trap axes at the Larmor frequency of 
r^L = 27r X 187 kHz. The field also produces a quadratic 
Zeeman energy shift qb = /i(70Hz/G^)5^ ^ h x 5 
Hz. The sample was also exposed to a linearly po- 
larized microwave- frequency magnetic field detuned by 
6 = ±27r X 40 kHz from the \F = l^mp = 0) to 
\F = 2^mF — 0) hyperfine transition. The resulting 
ac Zeeman shift of g^ = —Wt^j^^ was used to tune the 
total quadratic shift ^ = ^'b + ^/i by varying the Rabi 
frequency ^ and by choosing the sign of b [23l |24] . 

The spinor gas was then measured either by time-of- 
fiight analysis, through which we quantify the fractional 
Zeeman populations, ("m , of the "quasi-condensed" por- 
tion of the gas (defined conventionally by the bimodality 
of the spatial distribution of the expanded gas), or by 
in-situ measurements of the vector magnetization. We 
first discuss the time-of-fiight measurements, which were 
performed on gases that were prepared initially with frac- 
tional Zeeman populations (mj of either 0, 1/4 or 1/3. 

Evaporatively cooling such gases produces spinor con- 
densates containing around 3 x 10^ atoms at t = 0, and 
around 1.2 x 10^ atoms at t = 4 s. Assuming a con- 
ventional Thomas-Fermi distribution of condensed atoms 
within the optical trap, we determine the peak conden- 
sate densities to be no = 2.6 x 10^^ cm~^ at t = and 
1.8 xlO^^cm"^ at t = 4 s, respectively. We note that 
bimodal time-of-fiight distributions are observed already 
beginning at t c^ — 50 ms. At early times the spinor con- 
densates are found to have fractional Zeeman populations 
identical to those of non-degenerate gas from which the 
condensates form. Following that, spin- mixing collisions 
within the condensate cause the Zeeman populations to 
vary. For instance, starting with equal populations in 
the three Zeeman populations Fig. h[a), the ("ii popula- 
tion tends to decrease when the energy of the \mF = ±1) 
states is raised (^ > 0) , and to increase when it is lowered 
(g < 0). Comparing measurements on samples with dif- 
ferent initial spin compositions (Fig. [ijc)), we find that 
a common steady-state distribution among Zeeman lev- 
els is achieved within several seconds of equilibration, for 
\q\ ^ 10 Hz. 

We compare the measured Zeeman populations to the 
mean-field predictions for the uniform ground state of 
a ferromagnetic F = 1 spinor condensate [12 . For 
q < 0, the spin-dependent contact interaction is mini- 
mized by a fully magnetized condensate with magnetiza- 
tion aligned with the field axis x. For our experiments, 
the total magnetization (F^) = is conserved in the ab- 
sence of dipolar relaxation. To satisfy this constraint, 
one might expect the condensate to be divided into do- 
mains of magnetization oriented along ±x, resulting in 
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FIG. 1: Condensate Zeeman populations, Cii, during equili- 
bration of unmagnetized spinor gases at variable q. Measure- 
ments at various evolution times t = 200 ms (purple circles) , 1 
s (red squares), and 3 s (black diamonds) are shown with ini- 
tial C±i = 1/3 (a) and 1/4 (b). In (c), Zeeman populations at 

t = 2 s are compared for initial Cii = 1/3 (green triangles), 
1/4 (blue circles), and (yellow squares). For a narrow range 
M/h ^ 5 Hz (shaded), Zeeman populations evolve toward a 
common steady state within about 2 s, toward values qual- 
itatively in accord with mean-field predictions (solid line in 
c) based on the spin-dependent contact interaction. Outside 
that range, initial population differences persist throughout 
the experimentally accessed equilibration times. Data are av- 
erages of 3 experimental runs. 
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q < 2\c2\{n)y c^ h X 8 Hz, where 
{n)y = 1.85 X 10"^^cm~^ is the density averaged in the y 
direction after 2 s of equilibration, one expects that the 
condensate comprise transversely polarized spin domains 

(c) 

in which (^{ diminishes monotonically from 0.25 to zero 
[25] . At g = 0, fully magnetized domains of all spin 
orientations minimize the mean-field energy; one might 
therefore expect Cii = 1/3. These predictions ignore 
the number and composition of domain boundaries, and 
agree qualitatively with our measurements at small \q\. 

For larger values of 1^1, the Zeeman populations do not 
reach a common steady state within the 4 s lifetimes of 
our trapped samples. For q > 2|c2|no, the single- atom 
quadratic Zeeman shift is predicted to be sufficient to 
overcome the ferromagnetic spin-dependent contact in- 
teraction and lead to a condensate equilibrating purely 
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FIG. 2: Transverse (top) and longitudinal (bottom) magnetization after a variable evolution time at a quadratic shift of q/h = 
-5 (a), (b), and 5 (c) Hz for an initial spin mixture of Cii = 1/3. The transverse magnetization is represented using the 
color scheme shown, where the hue indicates the magnetization orientation and the brightness its magnitude. The maximum 
brightness is set by a a fully polarized condensate at 2 s evolution time. The longitudinal magnetization is represented by 
the color bar. Textures form from initially small domains and then coarsen to where only a few domains span the entire 
condensate. The late-time images show a spin-space anisotropy; for positive q the transverse magnetization appears brighter 
while for negative q the opposite trend is observed. 



in the \mF = 0) state p!2HT4] . However, the condensates 
prepared with ("ii 7^ retain non-equilibrium Zeeman 
populations nearly equal to their initial values, with the 
rate of Zeeman-population transfer diminishing with in- 
creasing q. For q < 0, the single-atom quadratic Zeeman 
shift favors population in the Imp = il) states. How- 
ever, for q/h < —20 Hz, while condensates within non- 
zero initial values of ("±1 do evolve toward the presumed 
equilibrium distribution, a condensate prepared in the 
non-equilibrium \mF = 0) Zeeman state shows almost 
no variation in its spin distribution. The extinction of 
spin-mixing dynamics at negative q for condensates in 
the \mF = 0) states was observed previously [23] . 

The strong variation of spin-mixing dynamics with 
variations in q that are far smaller than the thermal en- 
ergy (/cB(50nK) = h X 1 kHz) suggests that spin- mixing 
occurs primarily within the condensed fraction of the gas, 
with little influence of the non-condensed portion. The 
variation of condensate spin-mixing dynamics with q has 
been studied in the single-mode regime fT5l |26J , show- 
ing that the amplitude of Zeeman-population oscillations 
diminishes at large quadratic shifts. This reduction in 
amplitude occurs once the energy difference introduced 
between the two-atom states coupled by spin mixing be- 
comes larger than the spin-mixing interaction strength 
(l^l > 2|c2|n). However, it is unclear how this single- 
mode explanation applies to the inhomogeneous spin tex- 
tures studied here. 

While the Zeeman-population measurements are in- 
dicative of equilibration dynamics, they do not reveal 
directly the spontaneous magnetization of the quantum 
degenerate gas nor do they characterize the spatial distri- 
bution of the condensate spin. To explore these proper- 
ties, we employ magnetization-sensitive dispersive imag- 
ing. In contrast to previous work ^27l [28] , in which the 



magnetization was reconstructed from a series of images 
taken during Larmor precession, here we obtain a di- 
rect snapshot of the vector magnetization with just three 
images by applying spin-echo rf pulses between probes. 
Specifically, after a variable equilibration time, we apply 
the first imaging pulse, sending linear polarized probe 
light along the y axis and measuring its optical rota- 
tion to resolve one transverse component of the column- 
integrated magnetization M(p) with the position vector 
p in the imaged x-z plane. Following a delay tq of about 
50 /is, we apply a tt pulse with a resonant rf magnetic 
field. A second image, taken after an additional delay 
r = To + 7r/2QL, records the second component of the 
transverse magnetization. To measure the longitudinal 
magnetization component, we apply a properly timed tt- 
7r/2-7r pulse sequence before taking a third image. This 
spin-echo imaging sequence functions reliably even in the 
presence of > 27r x 10 kHz variations in Ql between ex- 
perimental runs. For future work, the relative parsimony 
of this method should allow either higher-resolution or 
repeated imaging of atomic spin distributions. 

The temporal evolution of spin textures produced upon 
cooling thermal mixtures with (^-^^ = 1/3 is shown in 
Fig.[2]for three different settings ofq. Within several hun- 
dred milliseconds, a highly corrugated spin texture de- 
velops within the condensate consisting of small domains 
with magnetization orientations distributed isotropically. 
The domains have a typical diameter of about 10 mi- 
crons, in agreement with previous observations [8 , al- 
though the textures in the present experiment do not 
show a preferred direction for spatial modulation of the 
magnetization in the x-z plane. 

The next few seconds of equilibration lead to distinct 
changes in both the composition and the spatial arrange- 
ment of the spin textures. The magnetization develops 



a g-dependent spin-space anisotropy visible in the con- 
trasting magnetization amplitudes for different quadratic 
shifts (Fig. [2]). Also evident is the development of larger 
domains of common magnetization. At early times the 
magnetization map consists of many randomly oriented 
spin domains with Neel-type domains walls, separating 
regions of common magnetization. Subsequently, within 
2-3 s, the texture contains just a few domains separated 
by Neel-type domains walls. 

To characterize these changes we consider the dimen- 
sionless correlation functions of the i component of the 
magnetization vector, given as 
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with /i being the atomic magnetic moment. This corre- 
lation function is evaluated over the central 35 x 85 /im^ 
area of the condensate. 

Setting Sp = 0, we obtain the variance of the longitudi- 
nal magnetization Gl{0)=Gx{0)^ and also the average of 
the transverse magnetization Gt(0) = {Gy{0)-\-Gz{0))/2 
[29] . For a fully magnetized gas Gl(0) + 2 x Gt(0) = 1. 
As shown in Fig. [3] (inset) for an initially unpolarized 
gas with an initial spin composition of C±i = 1/3, mag- 
netization at early equilibration times {t < 700 ms) de- 
velops with equal variance in the longitudinal and trans- 
verse orientations. Thereafter, within 2 s, the sum of 
the variances are consistent with unity and the magne- 
tization develops a spin-space anisotropy: easy-plane or 
easy-axis anisotropy seen from the enhanced transverse 
or longitudinal magnetization, for g' > or g' < 0, re- 
spectively (Fig. |3|. At q = 0^ the magnetization ori- 
entation remains isotropic. Both the g'-dependence and 
also the 1-2 s timescale for development of the spin-space 
anisotropy are consistent with the Zeeman-population 
measurements discussed above. 

A second trend in the evolution of these spin textures is 
the development of ever-larger domains of common mag- 
netization. The coarsening dynamics are characterized 
by measuring the area A of the central region of positive 
vector magnetization correlations, obtained by summing 
Gi{5p) over all spin components. The typical domain 
length / = y^AA/iT grows continuously from its initial 
value of about 10 jam to a final value / c:^ 40 /im over 
the 3 s equilibration time examined in our measurements 
(Fig. g. 

Phase-ordering kinetics from an initially disordered 
state have been studied in numerous systems [30 . Coars- 
ening dynamics is often found to yield a self-similar 
domain pattern with a characteristic length scale that 
evolves temporally as a power law l{t) ex t^, where the 
exponent b depends on the order parameter of the sys- 
tem. Adapting such classical theories to the coarsening of 
spin textures in a spinor condensate [311 [32] , one expects 
the conserved order parameter to vary effectively from a 
scalar order parameter (easy axis), for which b = 1/3, 
to a vector order parameter (easy plane/isotropic), for 
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FIG. 3: The magnetization variance versus q for C±i'^ = 1/3 
at t=2 s for transverse magnetization G'T(0)(red circles) and 
longitudinal magnetization Gl{0) (green triangles). For ^ = 
the longitudinal and transverse magnetization variances are 
equal, showing no preference in spin orientation (isotropic), 
while for ^ > (^ < 0), the data show a preference for trans- 
verse (longitudinal) magnetization, seen in the larger trans- 
verse (longitudinal) magnetization variance (easy plane/easy 
axis). (Inset) Temporal evolution of the transverse (G't(O), 
red circles) and longitudinal (G'l(O), green diamonds) mag- 
netization variance, evaluated within the central regions of 
spin textures produced from C±i — 1/3 spin mixtures at 
q/h = -5 Hz (left), Hz (middle) and 5 Hz (right). At 
early equilibration times, t < 700 ms, the magnetization has 
equal variance in the longitudinal and transverse orientations. 
Within 1-2 s, a clear preference for longitudinal orientation at 
^ > (middle), while for ^ = the spin distribution remains 
isotropic. Data are the average of 6 experimental repetitions, 
and error bars are the standard deviation of the data. 



which b = 1/4, according to the quadratic Zeeman shift 
j33l [34] . For g > we expect similar domain coarsening 
dynamics to the behavior observed for g = 0, whereas 
for g < 0, the rate of domain growth is expected to be 
slightly faster. 

The observed coarsening dynamics and the ability to 
tune the order parameter in a spinor condensate should 
enable a quantitative comparison with theories of phase- 
ordering kinetics. However, limitations specific to our 
present experiment prohibit a comparison with theory. 
These limitations include the finite size of the system, 
which when comparable to the domain size may modify 
the observed coarsening dynamics. 

In conclusion, we observe four temporal trends: a re- 
distribution of Zeeman populations, the development of 
magnetization from the initially unmagnetized gas, the 
development of a spin-space anisotropy, and the coars- 
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FIG. 4: Temporal evolution of domain length after variable 
evolution time at a quadratic shift of Hz for an (^±i = 1/3 
mixture. The domain length, as defined in the text, grows 
smoothly as a function of the evolution time, with domains 
of about 10 /im to final values of 40 /xm within 2 s. 



ening of spin domains. For small \q\^ all of these trends 
occur on a similar 1-2 s timescale, allowing the system 



to reach a common final state independent of the ini- 
tial spin composition. For larger |g|, equilibrium is not 
reached within the accessible experimental timescales. 
The timescale for spin mixing dynamics is set by the 
spin-dependent interaction energy, 2|c2|n ^ h x 8 Hz, 
resulting in spin mixing timescales of Tsm = h/{2\c2\n), 
far shorter than the observed seconds-long equilibration 
times. It remains unclear what determines the dynamical 
timescales for each of the observed processes. 
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